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Abstract 

We have used the stabihzed spin-polarized jeUium model to calculate the equilib- 
rium sizes of metal clusters. Our self-consistent calculations in the local spin-density 
approximation show that for an A'^-electron cluster, the equilibrium is achieved for 
a configuration in which the difference in the numbers of up-spin and down-spin 
electrons is zero or unity, depending on the total number of electrons. That is, a 
configuration in which the spins are maximally compensated. This maximum spin- 
compensation results in both the alternation in the average distance between the 
nearest neighbor ions and the odd-even alternations in the ionization energies of al- 
kali metal clusters, in a good agreement with the molecular dynamics findings and 
the experiment. These suggest a realistic and more accurate method for calculating 
the properties of metal clusters in the context of jellium model than previous jellium 
model methods. 
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1 Introduction 



The simplest model used in theoretical study of the properties of simple metal clusters 
is jellium model (JM) with spherical geometry. [|l], |^ In this model, the ions are 
replaced by a uniform positive charge background sphere of density n = 3/47rrf and 
radius R = (zNY^^Ts where z, N and are the valence of the atom, the number of 
constituent atoms of the cluster and the bulk value of the Wigner-Seitz (WS) radius of 
the metal, respectively. This model can be useful only when the pseudopotentials of 
the ions do not significantly affect the electronic structure. But, it is well-known that 
the JM has some drawbacks. |0, ^ Keeping the simplicity of the JM, the stabilized 
jellium model (SJM) of Perdew et a/. has lifted the essential deficiencies of the JM 
and significantly improved the results of calculations. In the SJM, the bulk metal (which 
is spin-unpolarized in the absence of external magnetic field) has been made stabilized 
through the introduction of a pseudopotential and fixing its core radius by a value that 
makes the pressure on the unpolarized bulk system to vanish. However, applying the SJM 
in the framework of rigid jellium background ^ may be suitable for closed-shell clusters 
of which the spin polarization of the valence electrons vanishes. It is well-known that, 
for example, the bond-length of a diatomic molecule depends on the relative orientations 
of the valence electrons. Hence, considering a metal cluster as a large molecule, or the 
bulk metal as a huge molecule, one should take account of the volume change due to spin 
polarization. We therefore, expect that the spherical jellium radius should be different for 
an A^-atom cluster with two different spin configurations. These facts led us to consider 
a stabilized jellium model in which the spin degrees of freedom be present. Stabilizing 
the jellium system with non-zero spin polarization, (, for the valence electrons resulted 
in the stabilized spin-polarized jellium model (SSPJM).]^ Here ( = (n^ — r;,j)/(n| + 
and rij, are the spin-up and spin-down electron densities, respectively. The SSPJM 
can be applied to metal clusters in two different ways. The first method, which has been 
used in Ref.Q, exploits the fact that the bulk metal expands as ( increases. We call that 
method as SSPJMl throughout this paper. In that method, the jellium sphere radius is 
taken as -R(C) = izNy^^fs{(). For fs{C) we had taken fs{C) = rs{0) + ^'^s(C) which 
fs(0) was the observed value of the bulk WS radius of the metal and Ar^ was obtained 
by the application of the local spin-density approximation (LSDA) to the infinite electron 
gas system. However, in our phenomenological accounting of the volume change, the 
core radius of the pseudo-potential for the electron-ion interaction has been considered 
as a parameter, which becomes polarization dependent as we force the pressure of the 
polarized bulk system to vanish. Using that scheme, we had calculated the energies of 
different metal clusters, both neutral and singly ionized, for different spin configurations 
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and had shown that instead of Hund's first rule for the ground state, the maximum spin- 
compensation (MSG) rule was governing. The MSG property which originates here from 
the polarization dependence of the core radius, leads to the odd-even alternations in the 
ionization energies that had been observed experimentally in the alkali metal clusters. 
On the other hand, if one assumes a fixed, polarization-independent form for the electron- 
ion interaction, the MSG property will be realized only for non-spherical geometries of 
the jellium background. 

We have recently shown p!0| that it is not always necessary for a finite spherical jellium 
system to increase its size as the polarization, (, is increased. This can be explained by 
considering the fact that for an open-shell cluster if one increases the spin polarization 
from the possible minimum value consistent with the Pauli exclusion principle, one should 
make a spin-flip in the last uncomplete shell. Because of high degeneracy for the spherical 
geometry, this spin-flip in the last shell does not change appreciably the kinetic energy 
contribution to the total energy but changes appreciably the exchange-correlation energy 
which in turn gives rise to a deeper effective potential that makes the Kohn-Sham (KS) 
1^ orbitals more localized and therefore a smaller size for the cluster. On the other 



hand, although the SSPJMl results in better ionization energies than the SJM [0] in that 
it reproduces the odd-even alternation, it always predicts incorrect cluster sizes. That 
is, in the SSPJMl, the equilibrium for a cluster is taken to be greater or equal to the 
bulk value of ( see Fig. 3 of Ref. |P] ), so that it approaches the bulk value from the 
above; whereas, the molecular dynamics (MD) results for the average distance between 
the nearest neighbor atoms show that the equilibrium value of the neutral clusters are 
less than the bulk value and it approaches the bulk value from the below [ see Fig 15(a) 
of Ref. |T^. To incorporate this correct behavior into our SSPJM calculations, which is 
the subject of this paper, we proceed parallel to the work of Perdew et a/.[jl2| for the spin- 
polarized case and call this method as SSPJM2. In the SSPJM2, for a given polarization, 
we first obtain the value of the core radius of the pseudopotential that stabilizes the bulk 
system, say r^, and then, using this value of in the energy functional of the cluster, 
we change the radius of the jellium sphere until the minimum energy is achieved. Our 
self-consistent calculations show that the absolute minimum energy corresponds to a spin 
configuration with maximum compensation as in the SSPJMl case. The equilibrium 
values corresponding to these minima lie below the bulk value, reproducing the correct 
behavior. These equilibrium values determine the equilibrium sizes of the clusters. If 
we plot the equilibrium value as a function of the number of constituent atoms in an 
alkali metal cluster, we see an alternating behavior, consistent with the MD results. [pT^ 
In this paper we have found the equilibrium properties of neutral and singly ionized Gs, 
Na and Al clusters of various sizes (2 < < 42) using jellium with sharp boundaries. We 
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have also repeated the SSP JM2 calculations using a jellium sphere with diffuse boundary. 



For the salce of comparison, we have derived the results of the work by Perdew et a/. [12 



which is denoted by SJMl. Comparing our SSPJM2 results with those of SSPJMl show 
that here, the average energies per electron and the ionization energies remain more or 
less the same but here, our SSPJM2 calculations show an improvement over the SSPJMl 
results for the equilibrium sizes of the clusters. 

In section |^ the calculational schemes has been explained. Section ^ is devoted to the 
results of our calculations and finally, we conclude this work in section |. 



2 Calculational Scheme 

In the context of the SSPJM, the average energy per valence electron in the bulk with 
density parameter and polarization ( is given by||^ 

€{rs, C) = ts{rs, C) + £xc{rs, C) + ^nirs, Tc) + euirs), (1) 

where ts and Sxc are noninteracting kinetic energy and exchange-correlation energies per 
electron, respectively, wr is the average value (over the WS cell) of the repulsive part of 



the Ashcroft empty core [ 14 1 pseudopotential. 



2z 2z 
w{r) = -—+WR, wr = +—0{r^-r), (2) 

and is given by wr = 3r^/rf. In Eq.(|^), z is the valence of the atom, and 9{x) is the 
ordinary step function which assumes the value of unity for positive arguments, and zero 
for negative values. The core radius, rc, will be fixed by setting the pressure of the bulk 
system equal to zero at equilibrium density n{() = 3/47rf^(C). In Eq.(|I]), is the average 
Madelung energy, = —9z/5ro. Here, ro is the radius of the WS sphere, tq = z^^^Vg, 
and for monovalent metals z = 1, and for polyvalent metals we set z* = 1 (for details 
see Ref.Q). All equations throughout this paper are expressed in Rydberg atomic units. 
The bulk stability is achieved when rc takes a value that makes the pressure of the system 
with a given ( to vanish at = fs{(): 



—€{rs,C,rc 
or,. 



= 0. (3) 

=rs{Q 

The derivative is taken at fixed ( and r^. Solution of the above equation gives the bulk 
value of Tc as a function of fg and (. Here, fsiC) is the equilibrium density parameter for 
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the bulk system with given ( and is evaluated by 

f,(C)=f,(0) + Arf«(C). (4) 

Here, fs{0) takes the observed value for a metal and for Cs, Na, and Al it takes the 
values of 5.63, 3.99, and 2.07, respectively. In the second term of the right hand side, 
the superscript "EG" refers to the electron gas, and Arf^i^O is evaluated by setting the 
pressure of the electron gas system equal to zero (see Eq.(19) of Ref.p]). The solution of 
Eq.(|^) at equilibrium density gives the bulk value of the core radius: 

^f(^.,C) = ^|-2t.(f.,C)-^x(r.,C)+r. (^J-^ ^c(r.,C)-^M(r.)| . (5) 

Now, using in the SSPJM energy functional of a cluster [ Eq.(20) of Ref.[§]], we obtain 
the SSPJM2 energy as 



+ {SvUs{r^,rs) I dv Q{v)[n{j) - n+(r)], (6) 



where 



+ ^y"c?r 0([n,n+];r)[n(r) - n+(r)] (7) 



and 



|,MiH-];r) = 2/* 



In Eq. (1), ((5f)ws is the average of the difference potential over the Wigner-Seitz cell 
and the difference potential, 5v, is defined as the difference between the pseudopotential 
of a lattice of ions and the electrostatic potential of the jellium background. 

The first and second terms in the right hand side of Eq.(|^) are the non-interacting 
kinetic energy and the exchange-correlation energy, and the last term is the Coulomb 
interaction energy of the system. In our spherical JM, we have 

n^{r) = ^^e{R-r) (9) 
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in which R = (zNY^^Ts is the radius of the jelhum sphere, and n(r) denotes the electron 
density at point r in space. 

Applying the SSPJM2 to an electron cluster with up-spin and down-spin 
electrons (A^ = A^^ + |) and polarization ( = {N^- — Ni)/{N-^ + N^^), the total energy 
becomes a function of A^, (, r^, and where is the density parameter of the jellium 
background and is given by Eq.(|^). The equilibrium density parameter, fs{N,(), for 
a cluster is the solution of 



0. 



=fs{N,0 

Here again, the derivative is taken at fixed values of A^, and r, 



(10) 



For an A^-electron 



cluster, we have solved the KS equations |]rT| self-consistently for various spin configura- 
tions and Vs values and obtained the absolute minimum-energy spin configuration and its 
corresponding density parameter. 



3 Results 

We have applied the SSPJM2 to calculate the equilibrium energies and sizes of different 
metal clusters. In our calculations for an A^ electron cluster, we have solved the KS 
equations for all possible spin configurations < C ^ 1 and obtained the minimum 
values of energies and corresponding values of each configuration. The self-consistent 
calculations for Cs, Na, and Al with 2 < N < 42 show that the absolute minimum- 
energy configuration obeys the MSG rule and the equilibrium value for the cluster, 
fs{N,(), is less than the bulk value because, for small clusters the ratio of surface to 
volume energies become comparable and the surface tension compresses the cluster. This 
effect is known as self-compression. In Fig. ^ we have compared the equilibrium 
values of "generic clusters", JMl (see Ref. [1^), with the SJMl results which reproduce 



correct trends. |12] To clarify the concept of the "generic cluster", suppose that one solves 
the self-consistent KS-LSDA equations for a spherical simple JM cluster with jellium 
radius equal to i? = N^^^Vg and total number of electrons A^. For a given A^, these 
calculations are performed for different values as well as different spin configurations 
until the equilibrium value, fs{N,(), corresponding to the absolute minimum-energy 



configuration is obtained ( See Fig. 4 of Ref.|10| ). Since in the calculations one does not 



use any specific parameter corresponding to a certain metal, the result does not simulate 
any real cluster, and we call it an A^-electron "generic cluster" . In the limit of A^ ^ oo, the 
infinite generic cluster tends to the electron gas system for which C and fg 4.18. As 
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is seen from the figure, the equihbrium Vg value for the generic cluster approaches the bulk 
value, 4.18, from the above which does not simulate the correct behavior for a real metal 
cluster. It is seen that for the generic clusters the equilibrium values are greater than the 
value rg = 4.18 whereas the SJMl predicts values that are smaller than the bulk value 
for Na (r^ = 3.99) in agreement with the MD findings. This comparison clearly shows 
that simple JM gives wrong molecular bond lengths. We have performed our calculations 
based on Eq. (|10|) both for jellium with sharp boundary, SSPJM2, and diffuse jellium, 
dif-SSPJM2. For our diffuse jellium calculations we have used the background density 
with the radial dependence as[|l5l 



( \ = l ^{1 - + t)e-^/'[smh{r/t)]/r}, r<R 
''+^''^ ~ \ n{l - ((i? + t)/2i?)(l-e-2«/*)}^e(^~'^/7r, r > R, ^ ' 

where n = 3/47rr^ , R = N^^^Vs, and t is a parameter related to the surface thickness. 
We have chosen t = 1 in all our diffuse jellium calculations and then have varied the 
value of Vg until the minimum energy is achieved. Figure 0(a) compares the equilibrium 
Vg values of neutral cesium clusters for different sizes. It is seen that in most cases (rather 
large clusters for which ^ << 1 ), the SSPJM2 and the SJMl predict the same values 
for the equilibrium r^; whereas for iV = 3, 5, and 7 the SSPJM2 predicts larger values. 
These larger values give rise to an alternation in the plot. The values obtained from the 
dif-SSPJM2 lie below the values obtained from the SSPJM2 and the SJMl. Figure |(b) 
shows the results obtained for singly ionized cesium clusters. In this case, we see that 
the alternations persist up to A^=15 and have relatively large amplitudes. For A^=3 i.e., 
singly ionized 4-atom cluster, the value obtained from the SSPJM2 has become larger 
than the bulk value which is related to the rough evaluation of Ar^. In Fig. §(a) we have 
shown the same quantities for neutral Na clusters. The behavior of SSPJM2 results is the 
same as in Fig. |^(a) but, in the case of the dif-SSPJM2 the value for A^=5 has become 
nearly equal to that of A^=6 and also, the value for A^=3 is less than that of N=4 which 
completely differs from the SSPJM2 results. Also, we could not find any finite value for 
N=2 case in the dif-SSPJM2. That is, as much as we decrease the input value, the 
total energy correspondingly decreases. This means that the surface tension dominates 
the internal pressure and collapses the cluster. Of course, this is not the case in reality and 
it is the consequence of the fact that here the surface thickness, t, has become comparable 
to the cluster radius, R. Figure ^(b) compares the plots of average distance between the 
nearest neighbors obtained from our SSPJM2 calculations and the MD calculations of 



Rothlisberger and Andreoni.|13] In order to estimate the average distance between the 



nearest neighbor ions in the cluster, we have assumed a bcc structure as in the bulk of 
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Na. Then the shortest distance between the ions, d, is related to the lattice constant, a, 
through d = a\/3/2. But, in the bcc structure for Na, there are two electrons in a cell 
and therefore, a = 2y^7r/3 r^. Combining these two relations results in d = y/n r^. The 
value Ts = 3.99 is appropriate for room temperature (T = 300 K) which results in a value 
of d = 7.07 bohrs whereas, for T = K the appropriate value for is 3.93 which gives 
rise to the value d = 6.96 bohrs. Therefore, our results should lie above the MD results 
[ see Fig 15(a) of Ref. |1^] because, the MD calculations have been performed for zero 
temperature. In Fig. ^(c) we have shown the plots of equilibrium Vg values for singly 
ionized sodium clusters as functions of number of electrons, N. The behavior is more or 
less the same as neutral one. Figure §(a) compares the SSPJM2 results for neutral Al 
clusters with the results obtained using the SJMl. Here, we have taken the effective value 
of z* = 1. The diamonds and squares in the plot show the physical points. The main 
difference between our results and the SJMl is in the size of the jellium atom of Al. In 
Fig. ^(b) we have compared the results for singly ionized Al clusters. The results show 
some differences for values of away from shell closings. Looking at the above-mentioned 
figures, we note that in all the three cases of Cs, Na, and Al the SSPJM2 results in a larger 
or equal values for the average distance than the SJMl, and in addition show alternations 
for small clusters. 

Finally, in Fig. ^(a) we have compared the plots of the total energies per electron in 
the two schemes of the SSPJM2 and the SSPJMl for Na clusters. It is seen that in the 
SSPJM2 the energies are relatively lower than those of the SSPJMl for smaller clusters 
but the same for larger ones. We have also calculated the ionization energies of the clusters 
using the dif-SSPJM2 and and compared with the dif-SSPJMl and experimental values 
in Fig. ^(b). Here also the odd-even alternations show up themselves in the SSPJM2 as 
well as in the dif-SSPJM2 results and the values obtained are more or less the same as in 
the dif-SSPJMl [see Fig. 7(c) of Ref. ||. Therefore, the SSPJM2 calculations for simple 
metal clusters has improved the previous work, SSPJMl, in that it not only reproduces 
the odd-even alternations in the ionization energies, but also it gives correct behavior for 
the equilibrium sizes of the clusters. 

4 Summary and Conclusion 

In this work we have performed the SSPJM calculations as in the case of ab initio molecular 
structure calculations. That is, we have firstly calculated the stabilizing core radius of 
the pseudopotential, rf, for the bulk system with nonzero spin polarization. Then, using 
this value in the energy functional of a cluster with given values of N and (, the energy 



8 



becomes a function of the single variable r^, the density parameter of the uniform jellium 
background. Minimizing this function with respect to gives us the equilibrium 
value and energy of the cluster with that specified N and (. Our self-consistent KS- 
LSDA calculations show that the equilibrium configuration is one in which the spins are 
maximally compensated as in our previous findings, This maximum spin compensation 
gives rise to the odd-even alternations seen in the experimental ionization energy plot of 
alkali metal clusters. Calculating the average distance between the nearest neighbors of 
Na clusters, we find a good agreement between our SSPJM2 results and those obtained 
from MD calculations. We have therefore improved our previous SSPJMl results in that 
the odd-even property is kept the same as before but here, the sizes of the smaller clusters 
have been predicted correctly. 
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Figure 1: The equilibrium values in atomic units as functions of the cluster size A^. 
The solid squares correspond to the "generic clusters" (JMl) defined in the text and the 
large squares correspond to Na clusters using the method of Ref.[T^. The dashed and solid 
lines correspond to the equilibrium rg value of the bulk (4.18) in simple JM and to the 
bulk value of sodium (3.99), respectively. 



Figure 2: (a) The equilibrium values in atomic units as functions of the cluster size for 
cesium clusters obtained from using the three schems of SJMl, SSPJM2, and dif-SSPJM2. 
The dashed line corresponds to the bulk value of for cesium (5.63). (b) Same as (a) 
but for singly ionized cesium clusters. 



Figure 3: (a) Same as Fig.^ for neutral Na clusters. The bulk value is 3.99. (b) The 
average distance between nearest neighbors in atomic units for Na clusters. The squares 
correspond to our findings through SSPJM2, appropriate for room temperature and the 
diamonds correspond to the molecular dynamics results at zero temperature. The dotted 
line correspond to the bulk value 7.07. (c) Same as (a) for singly ionized Na clusters. 



Figure 4: (a) The equilibrium values in atomic units for neutral Al clusters {z* = 1, 
Ts = 2.07). The diamonds and squares show the physical points in SSPJM2 and SJMl 
schemes, respectively, (b) Same as (a) for singly ionized Al clusters. 



Figure 5: (a) The total energies per atom of Na clusters in electron volts for the SSPJM2 
and the SSPJMl. The SSPJM2 results are somewhat lower than those of the SSPJMl 
for smaller clusters, (b) The ionization energies in electron volts for Na clusters in the 
dif-SSPJM2 and the dif-SSPJMl are compared with the experimental values. [^] 
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